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INTRODUCTION 

Several  researchers  have  produced  models  to  predict  the  flame 
spreading  and  pressurization  of  mobile  bed  of  solid  propellant.  A 
heated  debate  over  the  correctness  of  the  governlng^equatlons  has 
divided  Into  two  sides  - The  University  of  Illinois  model  against 
others  (Kuo^,  Gough^,  and  Fisher^).  In  addition  to  two  other  appli- 
cations'*^, Krler  and  associates^  have  applied  their  model  to  a 
BRL  laboratory  experiment®  which  detonated  a tetryl  pellet  at  one  end 
of  a circular  cylinder  filled  with  M30  propellant. 

It  Is  the  purpose  of  this  report  to  criticize  Krler 's  computa- 
tional results  and  to  present  the  different  results  predicted  by  Gough's 
model.  Direct  linkage  of  differences  In  theory  with  differences  In 
results  will  not  be  made  here  because  It  Is  still  too  complicated. 

This  report  will  provide  one  piece  of  limited  Information  on  the  two- 
phase  flow  modeling  debate. 


W,  VanTaesell  and  H.  Krier,  "Combustion  and  Flame  Spreading  Phenomena 
in  Gas-Permeable  Explosive  Materials",  Inti  J Heat  Mass  Transfer, 

18_  (12),  pl377-1386,  1973. 

K.  K.  Kuo,  J.  H.  Koo,  T.  R.  Davis,  and  G.  R.  Coates,  "Transient 
Combustion  in  Mobile  Gas  Permeable  Propellants" , Aata  Astronautiaa 
^ p573-S91,  1976. 

P.  S.  Gough,  "Numeriaal  Analysis  of  a Two  Phase  Flow  with  Explicit 
Internal  Boundaries",  Paul  Gough  Associates  PGA-TR-76-2,  September 
1976. 

E.  B.  Fisher,  K.  W.  Graves,  and  A.  P.  Trippe,  "Application  of  a 
Flame  Spread  Model  to  Design  Problems  in  the  155mm  Propelling  Charge", 
12th  JANNAF  Combustion  Meeting,  Newport,  RI,  CPIA  Publication  273, 
December  1975, 

M.  W.  Beckstead,  N.  L.  Peterson,  D.  T.  Pilcher,  B.  D.  Hopkins,  and 
H.  Krier,  "Convective  Combustion  Modeling  Applied  to  Deflagration-to- 
Detonation  Transition",  12th  JANNAF  Combustion  Meeting,  Newport,  RI, 
CPIA  Publication  273,  December  1975. 

H.  Krier,  S,  Rajan,  and  W.  F.  VanTassell,  "Flame  Spreading  and  Com- 
bustion in  Packed  Beds  of  Propellant  Grains",  AIAA  J 14J3),  March  1976. 

H.  Krier  and  S.  S.  Gokhale,  "Predictions  of  Vigorous  Ignition  Dynamics 
for  a Packed  Bed  of  Solid  Propellant  Grains",  Inti  J Heat  Mass 
Transfer,  19_  p91 5-923,  1976. 

J.  D.  Knapton,  I.  C.  Stobie,  and  R.  H.  Comer,  "Vigorous  Ignition  of 
M30  Propellant",  BRL  Memorandum  Report  No.  2662,  August  1976, 

AD  #B013322L. 
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It  has  recently  come  to  this  author's  attention  that  a similar  de- 
bate is  raging  in  the  two  phase  industrial  heat  transfer  community. 
Although  the  intensity  of  the  physical  processes  is  much  less  there, 
debate  on  formulation  of  the  conservation  equations  is  perhaps  more 
intense  than  in  the  gun  community.  For  example  see  Reference  9. 

EXPERIMENTAL 


Theoretical  Differences 


Conservation  equations  for  the  Illinois  model  are  derived  from 
the  principles  of  continuum  mechanics;  Gough's  equations  are  derived 
from  a formal  averaging  procedure  for  two  separate  phasesj^^  Term  by 
term  comparison  of  the  equations  has  been  reported  by  Kuo  . Although 
the  equations  for  mass  conservation  are  the  same  for  both,  the  equations 
for  momentum  and  energy  conservation  are  different. 


Controversial  "diffusion"  terms  in  Krier's  momentum  and  energy 
equations  arise  from  an  assumption  of  an  inviscid,  Newtonian  mixture. 
Having  assumed  inviscid,  Newtonian  components,  Gough  derives  no  such 
terms  in  his  equations.  Although  Beckstead  et  al^  showed  an  indiffer- 
ence to  the  presence  of  these  terms,  a detailed  examination  for  a gun 
problem  shows  these  terms  to  have  magnitudes  comparable  to  other  terms 
in  the  equations  during  at  least  part  of  the  ignition  sequence.  The 
equations  compare  as  follows: 

Gas  Momentum 


Krler:  - ^(^P)  - V + (u-u^)m^  + [(u-u^)^p(j]  - ij/u,  (1) 

Gough:  - 4)^  - V + (u-u  )m,  - i/<u  . (2) 

dx  P o 

Solid  Momentum 
Du 

Krier:  ^ ^ 

Gough:  -(u-u^)m^^  ^ ^ ~ ^ ~ (1-4')  ^ - (9) 


Proceedings  of  Two-Phaee  Flow  and  Heat  Transfer  Sympcsiwv-Workshop, 
October  1Q76,  Ft.  Lauderdale,  FL,  Clean  Energy  Research  Institute, 
University  of  Miami,  1976. 

K.  K.  Kuo,  "A  Summary  of  the  CANNAF  Workshop  on  Theoretical  Modeling 
and  Experimental  Measurements  of  the  Combustion  and  Fluid  Flow  Pro- 
cesses in  Cun  Prop>ellant  Charges",  13th  JA’JNAF  Combustion  Meeting, 
Monterey,  CA,  September  1976. 
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In  the  gas  phase  momentum  equation  the  difference  is  whether  the 
porosity  should  be  inside  or  outside  the  pressure  gradient.  In  the 
debate  about  the  form  of  the  pressure  gradient  term,  it  has  been  noted 
H 12 

by  Harlow  and  by  Gough  that  a term  would  accelerate  a system 

dx 

at  rest  with  a non-uniform  porosity.  Crowe  noted  that  in  a weight- 
less environment  where  such  a system  can  be  conceived,  there  should  be 
no  acceleration.  In  the  treatment  of  drag,  adding  Krler's  momentum 
equations  results  in  a net  drag  force  on  the  mixture  if  his  term  V 
is  the  same  in  both  equations.  Gough's  solid  phase  momentum  equation 
has  two  additional  forces  to  represent  stress  from  the  bed's  resistance 
to  compaction  and  a response  to  the  gas  pressure  gradient.  Having 
assumed  a fluidized  bed,  Krier  omits  the  stress  term. 

Similar  comparison  is  not  as  straight  forward  for  the  gas  energy 
equation  because  Krier  uses  total  energy  while  Gough  uses  internal 
energy  as  the  dependent  variable. 

DE  9 

Krier:  ^Pz—  = - (upi(i)  + 4ip  (u-u  )E  + (u-u  )4>p  - 

Dt  ox  mm 


3 2 

Ijip(u-u  ) + (fipu(u-u^)  - q^^  + 

m m gp 

(u-u 

^hem  - ^ -"-r^  ) • 

Gough's  form  can  be  converted  to  total  energy. 


(5) 


* DE 


+ ( 


2 2 
u -u 
2 


)] 


- 

^3x 


q + 


chem 


e + 


(6) 


There  are  several  differences  an  analysis  of  which  would  be  quite 
lengthy. 


F.  H.  Harlow  and  A.  A.  Amsden,  "Numerical  Calculations  of  Multi- 
phase Flow",  J Comp  Phys,  17_f  (1) , January  1973. 

P.  S.  Gough,  "The  Predictive  Capability  of  Models  of  Interior 
Ballistics",  12th  JANNAF  Combustion  Meeting,  Newport,  RI,  CPI A 
Publication  273,  December  1973. 

C.  T.  Crowe,  Washington  State  University , private  cormunication, 
1976. 
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No  such  parallel  comparison  arises  for  solid  phase  energy  conser- 
vation. Gough  (and  others^ »^)  treat  the  solid  only  as  a problem  in 
transient  heat  conduction  before  ignition.  When  the  computed  surface 
temperature  reaches  the  specified  ignition  temperature,  combustion 
takes  over  and  heat  is  exchanged  between  the  two  phases  only  by  the 
heat  release  of  combustion.  Krier  uses  an  equation  similar  to  that  of 
the  gas  to  conserve  energy  of  the  particles.  Solid  phase  bulk  tempera- 
ture is  the  dependent  variable  that  affects  heat  transfer  from  the  gas 
and  is  compared  to  some  criterion  for  Ignition.  In  a comparison  of  the 
surface  with  bulk  temperature  ignition  criterion,  reported  earlier^^,  a 
4K  bulk  temperature  rise  was  found  equivalent  to  a 155K  surface  tempera 
ture  rise  in  the  same  time  for  a gun  propelling  charge. 

The  Problem 


The  laboratory  experiment  loaded  a 25mm  diameter  chamber  190mm 
long  with  seven-perforated  M30  propellant.  After  a 38g  tetryl  pellet 
was  detonated  at  one  end,  pressure-time  was  measured  at  four  axial 
stations  in  the  bed  (63mm,  100mm,  138mm,  176mm). 

Code  input  values  are  given  in  Table  I.  Inputs  for  Gough's  code 
were  the  values  reported  by  Krier  except  where  noted. 

RESULTS  AND  DISCUSSION 


Pressure 


Figure  1 compares  the  predictions  of  pressure-time  at  two  bed 
locations  for  Krier 's  data  of  Reference  7 with  igniter  time  of  0.5ras. 

Comparison  of  these  predicted  results  with  the  experimental 
pressures  is  difficult  because  of  the  ragged  data.  In  six  tests,  peak 
pressures  ranged  from  40  to  300  MPa  at  the  63mm  gage,  and  130  to  220 
at  the  137mm  gage.  If  the  high  and  low  at  each  gage  were  disregarded, 
the  range  could  be  120  to  240  MPa  and  170  to  205  MPa  respectively. 

Estimates  for  the  speed  of  propagation  of  the  pressure  wave  in 
the  bed  vary  from  1 to  33mm/us  with  deceleration®.  Although  scatter 
of  the  data  precludes  firm  estimates,  the  predicted  rates  of  less  than 
Imm/ps  are  clearly  too  low.  The  experimental  rates  were  about  the 
same  for  inert  and  live  propellant  beds. 

Krier  reported  that  predicted  pressure  profiles  were  little 
changed  when  an  inert  propellant  was  assumed  (Fig.  8,  Ref.  7)  for  an 
igniter  time  of  0.1ms.  Gough's  code  shows  more  sensitivity  to  the 
difference.  The  dotted  lines  in  Figure  2 represent  inert  propellant 


C.  y.  Nelsori,  "Comparison  of  Predictions  of  Three  Two-Phase  Flow 
Codes",  ISth  JANNAF  Combustion  Meeting,  Monterey,  CA,  September  JIP76 


Table  I.  Code  Input  Data 


Propellant 


Outer  Diameter 

7.01mm 

Perforation  Diameter 

1.32mm 

Length 

15.85mm  ^ 

Density 

1.66g/cm 

Burning  Rate  Coefficient 

°-^°s(MPa)" 

Note  a 

Burning  Rate  Exponent 

0.667 

Weight 

99g 

Note  b 

Ignition  Temperature 
Speed  of  Sound  when 

450K 

Packed 

442m/s 

Settling  Porosity 

0.42 

Note  c 

Bed  Geometry 

Length 

200mm 

Initial  Temperature 

305K 

Initial  Porosity 

0.47 

Diameter 

25ram 

Thermodynamics 

Energy  Release  in 

Burning 

4736J/g 

Note  a 

Gas  Molecular  Weight 

23.2 

Gas  Specific  Heat  Ratio 

1.23 

Igniter  Energy 

3682J/g  3 

Note  d 

Covolume 

1.01  cm  /g 

Note  d 

Igniter  Discharge  Rate 

Rate 

17.84kg/cm/s 

Length 

12.5mm 

Notes : 

a.  Incorrectly  stated  in  Ref.  7. 

b.  Gough  Surface  Criterion  450K. 

c.  Not  used  in  Krier  code. 

d.  Not  stated  in  Ref.  7. 
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Figure  1.  Pressure  History  Comparisons,  Igniter  Time  0.5ms 


Pressure  History  Comparisons  for  Live  and  Inert  Propellant 


simulations.  Gough's  propagation  speed  for  the  pressure  wave  Is  about 
0.6  mm/us  while  Krler's  Is  about  0.3  mm/ys.  The  sound  speed  in  the  un- 
disturbed gas  is  about  0.3  mm/ys.  Test  data  reported  rates  1 to  3 mm/ys 
in  inert  beds.  Both  show  roughly  the  observed  proportional  drop  of 
about  half  in  peak  pressure  from  50  to  100mm  bed  position.  Although 
Krier  obtains  better  agreement  with  test  pressures  for  the  inert  bed, 
igniter  time  was  only  0.1  ms.  He  obtained  best  agreement  for  the  live 
propellant  test  results  with  0.5  ms  Igniter  time. 

Porosity 

From  the  theoretical  development,  Krier 's  porosity  profiles  should 
not  be  correct  wherever  the  bed  is  packed.  He  computes  porosity  in  the 
bed  center  as  0.250  when  the  initial  porosity  of  the  "packed"  bed  is 
0.470.  Such  a compression  cannot  be  co*"puted  with  a model  which  assumes 
the  bed  is  always  fluidized  with  no  particle  Interaction.  As  the  bed 
becomes  "more  packed"  the  propagation  of  disturbances  proceeds  through 
the  bed  as  though  it  were  a true  solid.  Propagation  rates  are  pro- 
bably Inversely  proportional  to  porosity.  Resistance  to  particle  motion 
increases  as  packing  increases  and  the  momentum  equation  must  account 
for  such  increased  resistance. 

The  Krier  code  sets  0.25  as  an  arbitrary  lower  limit  on  porosity. 
Packing  of  a bed  of  these  particular  grains  occurs  at  a porosity  above 
0.40.  For  any  porosities  below  the  0.40,  the  model  misrepresents  the 
mechanical  properties  of  the  bed.  Imposing  a lower  limit  has  the 
effect  of  creating  arbitrary  gradients  that  affect  the  coupled  equations. 
It  effectively  converts  solid  to  gas  without  combustion. 

Surprisingly,  the  porosity  predictions  by  the  two  codes  are  not 
far  apart  in  the  packed  bed  regions.  When  Gough's  packed  bed  sound 
speed  was  set  to  zero,  the  minimum  porosity  attained  was  about  0.25, 

Krier 's  lower  limit.  Predictions  of  other  flow  variables  by  Gough's 
code  were  little  changed  by  variation  in  packing  properties.  Pre- 
dictions of  each  code  are  shown  in  Figure  3. 

Gas  Temperature 

Computing  gas  temperatures  seems  a sensitive  task  in  such  numeri- 
cal calculations.  Numerical  problems  typically  show  themselves  in  os- 
cillations in  gas  temperature.  These  two  codes  predict  markedly  differ- 
ent gas  temperature  profiles  as  seen  in  Figure  4 for  a common  location 
of  the  ignition  front. 


Ih 


Figure  3.  Porosity  History  Comparisons 
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Figure  4.  Gas  Temperature  Comparisons 
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i-'fler's  temperature  profile  has  two  questionable  features  that 
need  ejcplanation.  At  the  igniter  end  temperature  seems  too  high.  Ig- 
niter gas  enters  at  a temperature  of  about  3400K  . Note 

that  energy,  not  enthalpy,  is  the  energy  added  in  the  equation.  With 
"N  no  compressive  work  on  the  gas  in  that  region,  the  temperature  should 

be  limited  by  the  temperature  of  the  entering  gas.  The  argument  is 
further  supported  by  later  results  in  the  same  calculation  which  pre- 
dict temperatures  over  lOOOOK.  A postulation  of  local  compression  is 
contradicted  by  the  monotonlc  increase  of  porosity  with  time.  If  the 
region  were  being  compressed,  as  is  the  region  ahead  of  the  pressuri- 
zation front,  local  porosity  would  decrease. 

At  the  front  of  the  compression  wave  in  the  bed  interior  Krler's 
predicted  gas  temperatures  and  heat  transfer  violate  thermodynamic 
principles.  The  initial  physical  condition  is  a quiescent  gas  in  ther- 
mal equilibrium  with  solid  particles.  Hot  gas  entering  at  the  aft 
portion  forms  a compression  front  driving  gas  and  particles  forward. 

In  the  forward  portion  of  the  bed  the  gas  should  be  heated  by  the  com- 
bination of  compression  by  particle  compaction  and  mixing  with  the  hot 
combustion  gas.  As  the  gas  temperatures  rise,  heat  is  transferred  by 
convection  (only  mode  allowed)  to  the  particles. 

What  the  code  predicts  however  is  a cooling  of  the  gas  from  305K  to 
about  140K  while  the  solid  phase  is  being  simultaneously  heated  from  305K 
to  310K.  Heat  transfer  from  a cold  gas  to  a hot  particle  is  inadmissable . 
Although  the  printed  output  in  the  appendix  of  Ref.  7 shows  only  one  time 
step,  the  full  results  show  the  minimum  temperature  region  propagating 
through  the  bed  but  never  any  particle  cooling.  Neither  the  low  gas  tempera- 
ture nor  the  particle  heating  can  be  justified  by  quantitative  arguments. 

Gough's  predictions  show  a more  credible  shape  of  temperature  dis- 
tribution although  the  magnitude  of  the  peak  seems  high.  There  is  a 
compression  wave  which  could  locally  heat  the  gas.  For  a temperature 
rise,  without  mass  addition,  to  1.4  times  the  uncompressed  temperature, 
tne  pressure  would  have  to  Increase  by  a factor  of  about  ten.  Local 
porosity  would  have  to  decrease  to  much  less  than  the  packing  porosity 
to  accomplish  such  a compression.  Computed  porosity  decrease  is  only 
a factor  of  two. 

When  an  excessively  large  igniter  source  was  assumed  (45  and  168 
kg/cm/s)  the  peak  temperatures  soared  to  values  far  in  excess  of  lOOOOK. 

Flame  Spreading 

Krler  reported  that  flame  spreading  rate  could  be  correlated  with 
assumed  igniter  source  strength  (Fig.  21,  Ref.  15).  For  the  mass 
addition  rate  used  in  this  comparison,  the  flame  spreading  rate  of 
0.3ram/iis  was  not  sensitive  to  the  Igniter  time.  Gough's  code  also 
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predicts  a dependence  of  flame  spreading  rate  on  Igniter  strength.  For 
three  values  used  (17.8,  143,  and  535  kg/cm/s)  the  corresponding  flame 
spreading  rates  are  (0.7,  1.1,  and  2.4  mm/ps) . Figure  5 shows  the  de- 
pendence for  both  Gough's  and  Krler's  codes.  The  rate  is  not  dependent 
on  the  time  duration  of  the  igniter.  Figure  6 compares  the  flame 
spreading  predictions  in  the  bed. 

The  dominant  Influence  of  igniter  strength  on  flame  spreading  rate 
can  te  explained  by  the  relative  gas  production  rates  of  igniter  and 
propellant.  Considering  the  grain  surface,  propellant  density,  and 
burning  rate,  combustion  gases  are  produced  at  a rate  of  0.041  P-67 
kg/cm/sec  which  at  100  MPa  is  only  0.9  kg/cm/sec  whereas  the  igniter 
input  is  17.8  kg/cm/sec.  As  the  wave  propagates  thru  the  bed  the  rate 
of  Increase  of  gas  in  the  region  behind  the  wave  is  controlled  by  the 
flow  of  gas  down  the  pressure  gradient.  In  a cell  where  the  pressure 
rise  rate  is  6 MPa/sec  (a  typical  computed  rate)  the  gas  increase 
rate  is  of  the  order  of  8 kg/cm/sec  of  which  only  about  0.9  kg/cm/sec 
is  provided  by  the  combustion.  Even  these  crude  estimates  show  the 
dominance  of  the  igniter  in  this  range  of  source  strengths. 

Phase  Velocities 

Krier  did  not  publish  plots  of  gas  and  solid  velocities.  Figures 
7,  8,  and  9 show  Gcjgh's  predicted  values  at  0.05,  0.11,  and  0.18  ms 
for  both  velocity  fields. 


CONCLUSIONS 

1.  Flow  predictions  are  different  for  the  Gough  and  Krier  codes. 

2.  Krier's  gas  temperature  predictions  in  the  bed  Interior  are  in- 
correct . 

3.  Both  codes  predict  a linear  dependence  of  flame  spreading  rate 
on  igniter  strength.  Gough's  rates  are  about  twice  Krier's. 


i 


20 


L 


V 


i 

I 


i 


IGNITER  STRENGTH  ( kg/cm/s  ) 


Figure  6.  Flainespreadlng 


TEMPERfllURE  OEO  K 


REFERENCES 


1.  W,  VanTassell  and  H.  Krier,  "Combustion  and  Flame  Spreading  Pheno- 
mena in  Gas-Permeable  Explosive  Materials",  Inti  J Heat  Mass 
Transfer,  (12),  pl377-1386,  1975. 

2.  K.  K.  Kuo,  J.  H.  Koo,  T.  R.  Davis,  and  G.  R.  Coates,  "Transient 
Combustion  in  Mobile  Gas  Permeable  Propellants",  Acta  Astronauti- 
ca  3 P573-591,  1976. 

3.  P.  S.  Gough,  "Numerical  Analysis  of  a Two  Phase  Flow  with  Explicit 
Internal  Boundaries",  Paul  Gough  Associates  PGA-TR-76-2,  September 
1976. 

4.  E.  B.  Fisher,  K.  W.  Graves  and  A.  P.  Trlppe,  "Application  of  a 
Flame  Spread  Model  to  Design  Problems  in  the  155mm  Propelling 
Charge",  12th  JANNAF  Combustion  Meeting,  Newport,  RI,  CPIA  Publi- 
cation 273,  December  1975. 

5.  M.  W.  Beckstead,  N.  L.  Peterson,  D.  T.  Pilcher,  B.  D.  Hopkins  and 
H.  Krier,  "Convective  Combustion  Modeling  Applied  to  Def lagratlon- 
to-Detonation  Transition",  12th  JANNAF  Combustion  Meeting,  Newport, 
RI,  CPIA  Publication  273,  December  1975. 

6.  H.  Krier,  S.  Rajan,  and  W.  F.  VanTassell,  "Flame  Spreading  and  Com- 
bustion in  Packed  Beds  of  Propellant  Grains",  AIAAJ  ^ (3)  p301-309, 
March  1976. 

7.  H.  Krier  and  S.  S.  Gokhale,  "Predictions  of  Vigorous  Ignition 
Dynamics  for  a Packed  Bed  of  Solid  Propellant  Grains",  Inti  J 
Heat  Mass  Transfer,  p915-923,  1976. 

8.  J.  D.  Knapton,  I.  C.  Stobie  and  R.  H.  Comer,  "Vigorous  Ignition  of 
M30  Propellant",  BRL  Memorandum  Report  No.  2662,  August  1976, 

AD  //B013322L. 

9.  Proceedings  of  Two-Phase  Flow  and  Heat  Transfer  Symposium-Workshop, 
October  1976,  Ft.  Lauderdale,  FL,  Clean  Energy  Research  Institute, 
University  of  Miami,  1976. 

10.  K.  K.  Kuo,  "A  Summary  of  the  JANNAF  Workshop  on  Theoretical 
Modeling  and  Experimental  Measurements  of  the  Combustion  and  Fluid 
Flow  Processes  in  Gun  Propellant  Charges",  13th  JANNAF  Combustion 
Meeting,  Monterey,  CA,  September  1976. 

11.  F.  H.  Harlow  and  A.  A.  Amsden,  "Numerical  Calculations  of  Multi- 
phase Flow",  J Comp  Phys  ^ (1),  January  1973. 


27 


I^EGEDI  -G  Page 


BL/.NK_vct  triL-KD 


REFERENCES 


12.  P.  S.  Gough,  "The  Predictive  Capability  of  Models  of  Interior 
Ballistics",  12th  JANNAF  Combustion  Meeting,  Newport,  RI , CPIA 
Publication  273,  December  1975. 

13.  C.  T.  Crowe,  Washington  State  University,  private  communication, 
1976. 


14.  C.  W.  Nelson,  "Comparison  of  Predictions  of  Three  Two-Phase  Flow 
Codes",  13th  JANNAF  Combustion  Meeting,  Monterey,  CA,  September 
1976. 

15.  H.  Krier  and  S.  Gokhale,  "Vigorous  Ignition  of  Granulated  Pro- 
pellant Beds  by  Blast  Impact",  BRL  Contract  Report  No.  263, 
October  1975,  AD  #B007026L. 


28 


DISTRIBUTION  LIST 


No . of 
Copies 


3 Organization  C 

Commander 

Defense  Documentation  Center 
ATTN:  DCC-TCA 

Cameron  Station 
Alexandria,  VA  22314 

Director 

Defense  Advanced  Research 
Projects  Agency 
ATTN:  C.R.  Lehner 

1400  Wilson  Boulevard 
Arlington,  VA  22209 

Director 

Institute  for  Defense  Analyses 
ATTN:  Dr.  H.  Wolfhard 

Mr.  R.T.  Oliver 
400  Army-Navy  Drive 
Arlington,  VA  22202 

Commander 

US  Army  Materiel  Development 
and  Readiness  Command 
ATTN:  DRCDMA-ST 

5001  Elsenhower  Avenue 
Alexandria,  VA  22333 


1 Commander 

US  Army  Aviation  Systems 
Command 

ATTN:  DRSAV-E 

12th  and  Spruce  Streets 

St.  Louis,  MO  63166 

1 Director 

US  Army  Air  Mobility  Research 
and  Development  Laboratory 
Ames  Research  Center 
Moffett  Field,  CA  94035 

1 Commander 

US  Army  Electronics  Connnand 

ATTN:  DRSEL-RD 

Ft.  Monmouth,  NJ  07703 


No . of 
Copies 


Organization 


Commander 

US  Army  Missile  Materiel 
Readiness  Command 
ATTN:  DRSMI-R 

Redstone  Arsenal,  AL  35809 

Commander 

US  Army  Tank  Automotive 
Development  Command 
ATTN : DRDTA-RWL 

Warren,  MI  48090 

Commander 

US  Army  Mobility  Equipment 

Research  & Development  Command 
ATTN:  Tech  Doc  Ctr,  Bldg.  315 

DRSME-RZT 

Ft.  Belvoir,  VA  22060 
Commander 

US  Army  Armament  Materiel 
Readiness  Command 
Rock  Island,  IL  61202 

Commander 

US  Army  Frankford  Arsenal 
ATTN:  DRDAR-LCE-CI , Dr.  J.  Lannon 

SCA-CC,  Mr.  C.  Dickey 
SCA-PP,  Mr.  L.  Stlefel 
Philadelphia,  PA  19137 

Commander 

US  Army  Armament  Research  and 
Development  Command 
ATTN:  DRDAR-LC,  Dr.  J.P.  Picard 

DRDAR-LCE,  Mr.  C.  Lenchitz 
DRDAR-LCE,  Dr.  R.F.  Walker 
Dover,  New  Jersey  07801 

Commander 

US  Army  White  Sands  Missile 
Range 

ATTN : STEWS-VT 

WSMR,  NM  88002 


DISTRIBUTION  LIST 


No.  of  No.  of 

Copies  Organization  Copies  Organization 


1 Commander  2 

US  Army  Watervllet  Arsenal 
ATTN:  R.  Thlerry/Code  SWEWV-RD 

Watervllet,  NY  12189 

1 Commander 

US  Army  Harry  Diamond  Labs  1 

ATTN:  DRXDO-TI 

2800  Powder  Mill  Road 
Adelphl,  MD  20783 

1 Commander  1 

US  Army  Materials  and  Mechanics 
Research  Center 
ATTN:  DRXMR-ATL 

Watertown,  MA  02172 

1 Commander 

US  Army  Natick  Research  and  2 

Development  Command 
ATTN:  DRXRE,  Dr,  D,  Slellng 

Natick,  MA  0.'762 

1 Director 

US  Army  TRADOC  Systems  1 

Analysis  Activity 
ATTN:  ATAA-SA 

WSMR,  NM  88002 


Commander 

US  Naval  Surface  Weapons  Center 
ATTN:  S.J.  Jacobs/Code  240 

Code  730 

Silver  Spring,  MD  20910 
Commander 

US  Naval  Surface  Weapons  Center 
ATTN:  Tech  Lib 

Dahlgren,  VA  22338 

Commander 

US  Naval  Underwater  Systems 
Center 

Energy  Conversion  Department 
ATTN:  R.S.  Lazar/Code  5B331 

Newport,  RI  02840 

Commander 

US  Naval  Weapons  Center 
ATTN:  Dr.  R.  Derr 

Mr.  C.  Thelen 
China  Lake,  CA  93555 

Commander 

US  Naval  Research  Laboratory 
ATTN:  Code  6180 

Washington,  DC  20375 


1 Commander 

US  Army  Research  Office 

ATTN:  Tech  Lib 

P.O.  Box  12211 

Research  Triangle  Park,  NC 

27706 


3 Superintendent 

US  Naval  Postgraduate  School 
ATTN:  Tech  Lib 

Dr.  David  Netzer 
Dr.  Allen  Fuhs 
Monterey,  CA  93940 


1 Chief  of  Naval  Research  2 

ATTN:  Code  473 

800  N.  Quincy  Street 
Arlington,  VA  22217 

1 Commander 

US  Naval  Sea  Systems  Command  2 
ATTN:  J.W.  Murrln  (NAVSEA-0331) 

National  Center,  Bldg.  2,  Rm  6E08 
Washington,  DC  20360 


Commander 

US  Naval  Ordnance  Station 
ATTN:  Dr.  A.  Roberts 

Tech  Lib 

Indian  Head,  MD  20640 
AFOSR 

ATTN:  J.F.  Mas! 

Dr.  B.T.  Wolfson 
Bolling  AFB,  DC  20332 


30 


DISTRIBUTION  LIST 


No.  of  No.  of 

Copies  Organization  Copies  Organization 


2 AFRPL  (DYSC) 

ATTN:  Dr.  D.  George 

Mr.  J.N.  Levine 
Edwards  AFB,  CA  93523 

1 Lockheed  Palo  Alto  Rsch  Labs 
ATTN:  Tech  Info  Ctr 

3521  Hanover  Street 
Palo  Alto,  CA  94304 

1 Aerojet  Solid  Propulsion  Co. 
ATTN:  Dr.  P.  Michell 

Sacramento,  CA  95813 

1 ARO  Incorporated 

ATTN:  Mr.  N.  Dougherty 

Arnold  AFS,  TN  37389 

1 Atlantic  Research  Corporation 
ATTN:  Dr.  M.K.  King 

5390  Cherokee  Avenue 
Alexandria,  VA  22314 

1 AVCO  Corporation 

AVCO  Everett  Research  Lab  Div 
ATTN:  D.  Stickler 

2385  Revere  Beach  Parkway 
Everett,  MA  02149 

2 Calspan  Corporation 

ATTN:  Dr.  E.B.  Fisher 

A.P.  Trlppe 
P.O.  Box  235 
Buffalo,  NY  14221 

1 ENKI  Corporation 

ATTN:  M.I.  Madison 

9015  Fulbrlght  Avenue 
Chatsworth,  CA  91311 

1 Foster  Miller  Associates,  Inc. 
ATTN:  A.J.  Erickson 

135  Second  Avenue 
Waltham,  MA  02154 


1 General  Electric  Company 
Armament  Department 
ATTN:  M.J.  Bulman 

Lakeside  Avenue 
Burlington,  VT  05402 

1 General  Electric  Company 

Flight  Propulsion  Division 
ATTN:  Tech  Lib 

Cincinnati,  OH  45215 

2 Hercules  Incorporated 

Alleghany  Ballistic  Lab 
ATTN:  Dr.  R.  Miller 

Tech  Lib 

Cumberland,  MD  21501 

2 Hercules  Incorporated 
Bacchus  Works 
ATTN:  Dr.  M.  Beckstead 

Dr.  R.  Simmons 
Magna,  UT  84044 

1 IITRI 

ATTN:  Dr.  M.J.  Klein 

10  West  35th  Street 
Chicago,  IL  60615 

1 Olln  Corporation 

Badger  Army  Ammunition  Plant 
ATTN:  J.  Ramnarace 

Baraboo,  WI  53913 

2 Olln  Corporation 

New  Haven  Plant 
ATTN:  R.L.  Cook 

D.W.  Riefler 
275  Winchester  Avenue 
New  Haven,  CT  06504 

1 Paiij.  Gough  Associates,  Inc. 
ATTN:  Dr.  P.S.  Gough 

P.O.  Box  1614 
Portsmouth,  NH  03801 


31 


DISTRIBUTION  LIST 


No . of 

Copies  Organization 

1  Physics  International  Company 
ATTN:  J,D.  Watson 

2700  Merced  Street 
Leandro,  CA  94577 

1 Pulsepower  Systems,  Inc. 

ATTN:  L.C.  Elmore 

815  American  Street 
San  Carlos,  CA  94070 

1 Science  Applications,  Inc. 

ATTN:  Dr.  R.B.  Edelman 

23146  Cumorah  Crest 
Woodland  Hills,  CA  91364 

2 Rockwell  International  Corp. 
Rocketdyne  Division 

ATTN:  Dr.  C.  Obert 

Dr.  J.E.  Flanagan 
6633  Canoga  Avenue 
Canoga  Park,  CA  91304 

2 Rockwell  International  Corp. 
Rocketdyne  Division 

ATTN:  Mr.  W.  Haymes 

Tech  Lib 

McCJregor,  TX  76657 

1 Shock  Hydrodynamics,  Inc. 
ATTN:  Dr.  W.H.  Anderson 

4710-16  Vineland  Avenue 
North  Hollywood,  CA  91602 

1 Thlokol  Corporation 
Elkton  Division 
ATTN:  E.  Sutton 

Elkton,  MD  21921 

3 Thiokol  Corporation 

Huntsville  Division 
ATTN:  Dr.  D.  Flanigan 

Dr.  R.  Click 
Tech  Lib 

Huntsville,  AL  35807 


No . of 

Copies  Organization 

2 Thlokol  Corporation 
Wasatch  Division 
ATTN:  Dr.  John  Peterson 

Tech  Lib 
P.O.  Box  524 

Brigham  City,  UT  84302 

1 TRW  Systems  Group 

ATTN:  Dr.  H.  Korman 

One  Space  Park 

Redondo  Beach,  CA  90278 

2 United  Technology  Center 

ATTN:  Dr.  R.  Brown 

Tech  Lib 
P.O.  Box  358 
Sunnyvale,  CA  94088 

1 Universal  Propulsion  Co. 
ATTN:  H.J.  McSpadden 

P.O.  Box  546 
Riverside,  CA  92502 

1 Battelle  Memorial  Institute 
ATTN:  Tech  Lib 

505  King  Avenue 
Columbus,  OH  43201 

1 Brigham  Young  University 

Dept  of  Chemical  Engineering 
ATTN:  Prof,  R.  Coates 

Provo,  UT  84601 

1 California  Institute  of  Tech 
Jet  Propulsion  Laboratory 
ATTN:  Prof,  F.E.C.  Cullck 

Pasadena,  CA  91103 


32 


/ 


DISTRIBUTION  LIST 


No.  of 

Copl'ea  Organisation 

1 Case  Western  Reserve  Univ. 

Division  of  Aerospace  Sciences 
ATTN:  Prof.  J.  Tien 

Cleveland,  OH  44135 

3 Georgia  Institute  of  Technology 
School  of  Aerospace  Engineering 
ATTN:  Prof.  B.T.  Zlnn 

Prof.  E.  Price 
Prof.  W.C.  Strahle 
Atlanta,  GA  30332 

1 Institute  of  Gas  Technology 
ATTN:  Dr.  D.  Gidaspov 

3424  S.  State  Street 
Chicago,  IL  60616 

1 Johns  Hopkins  University 
Applied  Physics  Laboratory 
Chemical  Propulsion  Informa- 
tion Agency 

ATTN:  Mr.  T.  Christian 

Johns  Hopkins  Road 
Laurel,  MD  20810 

1 Massachusetts  Institute  of 
Technology 
Dept  of  Mechanical  Engineering  1 
ATTN:  Prof.  T.  Toong 

Cambridge,  MA  02139 

1 Pennsylvania  State  University 
Applied  Research  Lab 
ATTN:  Dr.  G.M.  Faeth 

P.O.  Box  30 

State  College,  PA  16801 

1 Pennsylvania  State  University 

Dept  of  Mechanical  Engineering  1 
ATTN:  Prof.  K.  Kuo 

University  Park,  PA  16801 


No.  of 

Copies  Organization 

3 Princeton  University 

Forrestal  Campus  Library 
ATTN:  Dr.  L.  Caveny 

Tech  Lib 

Prof  M.  Summerfield 
P.  0.  Box  710 
Princeton,  NJ  08540 

2 Purdue  University 

School  of  Mechanical  Engineering 
ATTN:  Prof.  J.  Osborn 

Prof.  S.N.B.  Murthy 
TSPC  Chaffee  Hall 
West  Lafayette,  IN  47906 

1 Rutgers  State  University 
Dept  of  Mechanical  and 
Aerospace  Engineering 
ATTN:  Prof.  S.  Temkln 

University  Heights  Campus 
New  Brunswick,  NJ  08903 

1 Southwest  Research  Institute 
Fire  Research  Section 
ATTN:  W.H.  McLain 

P.O.  Drawer  28510 
San  Antonio,  TX  78228 

Stanford  Research  Institute 
Propulsion  Sciences  Division 
ATTN:  Tech  Lib 

333  Ravens wood  Avenue 
Menlo  Park,  CA  94024 

1 Stevens  Institute  of  Technology 
Davidson  Laboratory 
ATTN:  Prof.  R.  McAlevy,  III 

Hoboken,  NJ  07030 

University  of  California, 

San  Diego 
AMES  Department 
ATTN:  Prof.  F.  Williams 

P.O.  Box  109 
La  Jolla,  CA  92037 


» 


.J 


33 


DISTRIBUTION  LIST 


No . of 

Copies  Organization 

1 University  of  Illinois 

Dept  of  Aeronautical  Engineering 
ATTN:  Prof.  H.  Krler 

Transportation  Bldg,  Rm  105 
Urbana,  IL  61S01 

1 University  of  Minnesota 

Dept  of  Mechanical  Engineering 
ATTN:  Prof.  E.  Fletcher 

Minneapolis,  MN  55455 

2 University  of  Utah 

Dept  of  Chemical  Engineering 
ATTN:  Prof.  A.  Baer 

Prof.  G.  Flandro 
Salt  Lake  City,  UT  84112 

1 Washington  State  University 

Dept  of  Mechanical  Engineering 
ATTN:  Prof.  C.  T.  Crowe 

Pullman,  WA  99163 

Aberdeen  Proving  Ground 


Marine  Corps  Ln  Ofc 
Dir,  USAMSAA 


